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Abstract 

The roots of Bernoulli polynomials, B n (z), when plotted in the complex plane, accumulate around 
a peculiar H-shaped curve. Karl Dilcher proved in 1987 that, on compact subsets of C, the Bernoulli 
polynomials asymptotically behave like sine or cosine. Here we establish the asmptotic behavior of 
B n (nz), compute the distribution of real roots of Bernoulli polynomials and show that, properly rescaled, 
the complex roots lie on the curve e _2?rIm ' z ' = 2ire\z\ or e 2 ^ lm ^ = 2ne\z\. 



This paper came out of the author's Geometry project at Perm State's Mathematics Advanced Study 
Semesters, Fall 2006. Adrian Ocneanu showed his geometry class pictures suggesting the roots of Bernoulli 
polynomials lie on a distinct curve. See figure 1. To prove this result, the author learned about the asymptotic 
behavior of the Bernoulli polynomials. Observing convergence behavior by plotting the roots of a sequence 
of polynomials is not new. Gabor Szego showed in 1928, the roots of P n (z) — YHk=o fj ^ e on ^ ne curve 
e ReO) _ e | z |^ j ean Dieudonne read Szego's paper and prove it a different way in 1935. Also there is a AMS 
memoir by Karl Dilcher from 1987, Dil], where he looks for parabolic zero-free regions in the complex plane. 
In this paper we find the exact curve around which the rescaled roots of Bernoulli polynomials cluster. 



1 The Bernoulli Polynomials 

Theorem 1.1. The roots of B n (nz) lie on the curve e ~ 2 ^ Im i z ) — 2ire\z\ i/Im(z) > and e 27rIrn ( z ) = 2ire\z\ 
iflm(z) < 0. 



Figure 1: The roots of Biqq(z) — via Mathematica 



Let's quickly review some definitions and properties of the Bernoulli numbers and their polynomials. The 
Bernoulli numbers are defined by a generating function: 
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There is no simple way to generate the Bernoulli numbers, they are defined as the coefficients of the Taylor 
expansion of the function above. The Bernoulli numbers are also coefficients of the Bernoulli polynomials 
which are also defined by a generating function in two variables: 
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One use of the Bernoulli polynomials is to express sums of similar powers: 

N 

{N+l)Y / k m = B m+1 (N + l)-B m+1 (1) 

k=0 

In fact Ocneanu's original question as posed in class was to consider the left side of (1), consider it as a 
polynomial in a complex variable and plot its roots. This his question was about the sums of like power 
and not necessarily Bernoulli polynomials. However, we simplify the problem silghly by not subtracting off 
terms. 



2 History of the Problem 

Given a sequence of polynomials, how do we know the roots accumulate on any curve? Futhermore, how 
do we compute this curve? In the October 2005 edition the American Mathemaitcal Monthly I saw this 
plot: 




Figure 2: The roots of Y^k=a ^k^ = fit the curve similar to \ze z x | = 1 



The paper, |Zem| . was by Stephen Zemyan of Pennsylvania State University, Mont Alto. He wrote on the 
roots of the truncated exponential series. If you cut the taylor series of e x after N terms you get a polynomial 

N-i k 
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As N — ► oo the roots go to oo since e x has no complex roots, but if you rescale by a factor of N, the roots 
of Pn(Nx) approach the curve |ze z_1 | = 1. This was proven by Gabor Szego in 1928 appearing in his 
famous problem book |Sze2j and in a seperate paper [Szej . Basically Szego's result has to do with the gamma 
function that 

k=0 u 

However, this work seemed to come out of nowhere and I went in search of more systemattic approach to 
this problem. 
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Around 1935, Jean Dieudonne proved the same result about the rescaled roots of the truncated exponential 
series. He shows that: 

n\[e nz - P n (nz)} 1 



[1 + Xn(z)} 



(2) 
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where X n (z) — > uniformly on the {z : \z\ < 1, \z — 1| > e}. Bounding away from 1 is essential because 
of the pole in the asymptotics. In his paper Dieudonne says that his solution is a technique which works 
for general classes of functions. Having read Dieudonne's proof, the author will try this technique on the 
Bernoulli polynomials. 



Unlike the exponential problem proved by Szego and Dieudonne, the Bernoulli polynomials are not the 
truncations of any particular infinite series. The coefhents seem to "jump" around and diverge. There was no 
obvious way to get around this problem In 1987, Karl Dilcher showed the following formulas using properties 
of the Ricmann Zeta funtion: 

Theorem 2.1. ,4s in JDWj . let 



T 2k (z) = Y(-l) 
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idT 2k+1 ( Z ) = j2(-iy 
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be the truncated sine and cosine series, for n € N. Then 



C-ir<g£«.(. + i)-r..,:. 



< 2- n e 47r|z| 



Corollary 2.1. On compact subsets of C we get uniform convergence: 

cos (27rz) if n is even 

sin (2ttz) if n is odd 



lim ^B n {z) 



This showed the Bernoulli polynomials were intimately related to the basic trigonometric functions and 
gave a quantative error between the limits of the respective truncations. However, if we substitute nz as an 
argument instead of z, we get uniform convergence on a small disc where e 47r ' z '/2 < 1 or \z\ < ln2/47r — 
0.055 < 0.159 = 1/27T. This means we cannot show these functions are unformly close on the full disc 
on which the roots lie. However, we do get a compelling picture 
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Figure 3: The roots of the Bernoulli polynomials compared to the locus of roots of truncated cosine. 



3 The Real Case 



Now we will show the real roots of B n (nz) are uniformly spaced on an interval [— l/27rne, \j2-nne\. Since 
Karl Dilcher related Bernoulli polynomials to trigonometric series and exponential functions, Fourier series 
may be the appropriate tool. Let's recall a result of Hurwitz: 
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Lemma 3.1. The Fourier series of the Bernoulli polynomials are given by: 



v ; fe >0 



Notice this Fourier series has period 1 so we use the fractional part sign {x}. 
Proof. The coefficients can be found all at once using the generating function: 
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For each individual fe we get a Fourier series like the one above. □ 

What are the asymptotics of this periodized version of the Bernoulli polynomials? 

Lemma 3.2. The Bernoulli polynomials B n ({x}) converge uniformly to 

(2ni) n e 27rix + e- 2lTlx 
x 



2n\ 2 

in the unit disc D. 

Proof. Consider the Fourier series from the previous lemma: 

n! e 2irikx 

y 1 |fe|>o 

Every term except k — ±1 goes to zero as n goes to infinity. This gives us pointwise convergence. Is this 
convergence uniform for x £ D? Yes 
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Now we have estimated the periodized version of the Bernoulli polynomials B n ({z}) but how do we estimate 
B n (mx) - B n ({mx})l 

Lemma 3.3. As n approaches infinity B n {mx) — B n ({mx}) approaches 

nm 71 ' 1 f mt n_1 dt 
Jo 

This convergence is uniform in the unit interval [0, 1 — e] for any e > 0. 
Proof. Let's do one at a time: B n (x + 1) — B n {x) = nx 11 ^ 1 . Then 

\nx\ — 1 

B n (mx) - B n ({mx}) = ^ B n ({mx} + k + 1) - B n ({mx} + k} 
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\_nx\ — 1 
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= nm x — > m 
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nm 11 - 1 x / mt n - x dt 
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For x G [0, 1], does the Riemann sum above convergence uniformly in B to the following integral? 

mt n - l dt = x n 



Since x m is monotone in x on [0, 1], the difference between the upper and lower Riemann sums is just the 
difference between the first term of the upper sum and the last term of the lower sum. The Riemann sum 
above is the lower sum. The upper sum is defined between k — 1 and k = \nx\ . Their difference is just 
x™" 1 — [{mx} /m] n ~ 1 < x™ -1 . This means convergence will be uniform for x G [0, 1 — e] for e > 0. □ 



How let's take another look at (3): 
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We need an estimate of Stirling's formula with errors. By considering the integral 
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e t(ln x -x) da 



it is possible using Laplace's method to get Stirling's esimate 
Thus 
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Asymptotically the roots of B n (nx) — B n ({nx}) + (B n (nx) — B n ({nx})) are the roots of e 2 
This is basically the stategy we will use in passing to the complex roots. 



4 The Complex Case 

The previous section's result can be found in [Vcs from 1999. In this section, while the proof is new, the 
result was anticipated by Karl Dilcher. Let's start writing the simple generating function relation 

f xt f (x+yi)t 

Jl_ x e yit = 



and rewriting it in terms of Bernoulli polynomials: 

n 

B n (x + yi) = J2(l) B k(x)(yi) 



fc=0 

Let's define a functional: 



n — k 



oo n co „ " / \ 

5>„(.x)- - e^]TF„(*)- or F n (x) -> £ ("WsXyi)"-* 

k=0 ' k=0 ' fc=0 ^ ' 



This functional is linear in generating functions F{x,t) or equivalently, sequences F n . In particular, we 
evaluate this functional for F n (x) = B n ({nx}) and F n {x) = B n (nx) — B n {{nx\). 
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Lemma 4.1. As n approaches infinity, B n (nx + niy) — B n ({nx} + iny) approaches 

n n [ nix + yiy^dx 
Jo 

Proof. 



J2(l)(B k (nx)-B k ({nx m y^ = „» ± Q x \ ^ \ (^^) ^ 
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This cooresponds to a Riemann integral: 



n(x + yi)" da; 

and since this function is monotone in x we get the same error bound as before: 

n(x + yi)™- 1 - n f + y j 



As long as \x + yi\ < 1 — e wc get uniform convergence for e > 0. 

Now let's give our final major estimates of the paper: 
Lemma 4.2. As n approaches infinity B n ({nx} + iny) approaches 

[ n O( n l/2"(l + !/„))] 



□ 



(27rie 



[P n (2imy) + P n (-2nny)} 



Proof. Now we transform B n ({nx}). We already know its Fourier series expansion for real arguments. Let's 
use the convolution formula above to get add an imaginary value iny : 



B n ({nx} + iny) = ^ I " j B m ({ra})(my) n " 

m=0 
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Therefore we get the following estimate for B n (nx + iny). We pull out the |fc| = 1 terms in the Fourier series, 
which do not decay with m. In the third line, we use our uniform bound for the |fc| > 2 tail in Lemma 3.2 
to give line four. 



B n ({nx} + iny) 



in-m)\ ^ 
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Having summed all the errors we get a worrisome error which might grow as ■^0(P n {±Aimy)) which might 
grow faster than the dominant terms! Indeed our proof would collapse if we estimated this as ±0(e ±47rny ), 
but it so happens that ^ P n {Aim\y\) < P n (2wn\y\). Let's double-check this: 



1 1 n_1 

Pn( 4nn\y\) = _ £ 



(47m|y|) A 
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= P n (27rn\y\) 
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Now we use Dieudonne's estimate to say that within the unit circle our Taylor series converges to 
P n (2irinz) — > e 2mnz uniformly inside the disc B/27T. 

n\[e 2mnz ~ P„(2*in*)] = [1 + X n (z)} 

This is to say the P n (2-Kinz) converges to the exponential uniformly on the disc j^B. Also note that 
[x n {l + O(l)] 1 /™ = x[l + 0(l/n)}. 



B n ({nx} + iny) 
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and we use Stirling's formula to get the above estimate. 
Earlier we showed: 

B n (nx + niy) — B n ({nx} + iny) = [n(x + iy)] r ' 

This tells us this equation we are looking for is: 

) . _ , e 2 ™y if y > 
""" • ~ ' e~ 2 ™y if y < 



x + iy 
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